ACCRETION DISC THEORY: FROM THE 
STANDARD MODEL UNTIL ADVECTION 



G.S.BISNOVATYI-KOGAN* 



Abstract 

Accretion disc theory was first developed as a theory with the local 
heat balance, where the whole energy produced by a viscous heating was 
emitted to the sides of the disc. One of the most important new inven- 
tion of this theory was a phenomenological treatment of the turbulent 
viscosity, known as "alpha" prescription, when the (r<j>) component of the 
stress tensor was approximated by (a P) with a unknown constant a. 
This prescription played the role in the accretion disc theory as well im- 
portant as the mixing-length theory of convection for stellar evolution. 
Sources of turbulence in the accretion disc are discussed, including non- 
linear hydrodynamical turbulence, convection and magnetic field role. In 
parallel to the optically thick geometrically thin accretion disc models, 
a new branch of the optically thin accretion disc models was discovered, 
with a larger thickness for the same total luminosity. The choice between 
these solutions should be done of the base of a stability analysis. The 
ideas underlying the necessity to include advection into the accretion disc 
theory are presented and first models with advection are reviewed. The 
present status of the solution for a low-luminous optically thin accretion 
disc model with advection is discussed and the limits for an advection 
dominated accretion flows (ADAF) imposed by the presence of magnetic 
field are analysed. 

1 Introduction 

Accretion is served as a source of energy in many astrophysical objects, includ- 
ing different types of binary stars, binary X-ray sources, most probably quasars 
and active galactic nuclei (AGN). While first developement of accretion theory 
started long time ago (Bondi and Hoyle, 1944; Bondi, 1952), the intensive devel- 
opement of this theory began after discovery of first X ray sources (Giacconi et 
al, 1962) and quasars (Schmidt, 1963). Accretion into stars, including neutron 
stars, is ended by a collision with an inner boundary, which may be a stellar 
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surface, or outer boundary of a magnctosphcrc for strongly magnetized stars. 
We may be sure in this case, that all gravitational energy of the falling matter 
will be transformed into heat and radiated outward. 

Situation is quite different for sources containing black holes, which are dis- 
covered in some binary X-ray sources in the galaxy, as well as in many AGN. 
Here matter is falling to the horizon, from where no radiation arrives, so all 
luminocity is formed on the way to it. The efficiency of accretion is not known 
from the beginning, contrary to the accretion into a star, and depends strongly 
on such factors, like angular momentum of the falling matter, and magnetic 
field embedded into it. It was first shown by Schwarzman (1971), that during 
spherical accretion of nonmagnetized gas the efficiency may be as small as 10~ 8 
for sufficiently low mass fluxes. He had shown that presence of magnetic field 
in the accretion flux matter increase the efficiency up to about 10%, and ac- 
count of heating of matter due to magnetic field annihilation in the flux rises 
the efficiency up to about 30% (Bisnovatyi-Kogan, Ruzmaikin, 1974). In the 
case of a thin disc accretion when matter has large angular momentum, the ef- 
ficiency is about 1/2 of the efficiency of accretion into a star with a radius equal 
to the radius of the last stable orbit. Matter cannot emit all the gravitational 
energy, part of which is absorbed by the black hole. In the case of geometrically 
thick and optically thin accretion discs the situation is approaching the case of 
spherical symmetry, and a presence of a magnetic field playes also a critical role. 

Here we consider a developement of the theory of a disk accretion, start- 
ing from creation of a so called "standard model", and discuss recent trends, 
connected with a presence of advection. 

2 Developement of the standard model of the 
disc accretion into a black hole 

Matter falling into a black hole is gathered into a disc when its angular momen- 
tum is sufficiently high. It happens when the matter falling into a black hole 
comes from the neighbouring ordinary star companion in the binary, or when 
the matter appears as a result of a tidal disruption of the star which trajectory of 
motion approaches sufficiently close to the black hole, so that forces of selfgrav- 
ity could be overcomed. The first situation is observed in many galactic X-ray 
sources containing a stellar mass black hole (Cherepashchuk, 1996). A tidal 
disruption happens in quasars and active galactic nuclei (AGN), if the model 
of supermassive black hole surrounded by a dense stellar cluster of Lynden-Bell 
(1969) is true for these objects. 

The models of the accretion disc structure around a black hole had been in- 
vestigated by Lynden-Bell (1969), Pringle and Rees (1972). The modern "stan- 
dard" theory of the disc accretion was formulated in the papers of Shakura 
(1972), Novikov and Thorne (1973) and Shakura and Sunyaev (1973). It is im- 
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portant to note, that all authors of the accretion disc theory from USSR were 
students (N.I.Shakura) or collaborators (I.D.Novikov and R.A.Sunyaev) of aca- 
demician Ya.B.Zeldovich, who was not among the authors, but whose influence 
on them hardly could be overestimated. 

The equations of the standard disc accretion theory were first formulated by 
Shakura (1972); some corrections and generalization to general relativity (GR) 
were done by Novikov and Thorne (1973), see also correction to their equations 
in GR made by Riffert & Herold (1995). The main idea of this theory is to 
describe a geometrically thin non-self-gravitating disc of the mass Ma, which is 
much smaller then the mass of the black hole M, by hydrodynamic equations 
averaged over the disc thickness 2h. 



2.1 Equilibrium equations 

The small thickness of the disc in comparison with its radius h <C r indicate to 
small importance of the pressure gradient VP in comparison with gravity and 
inertia forces. That leads to a simple radial equilibrium equation denoting the 
balance between the last two forces occuring when the angular velocity of the 
disc Q is equal to the Keplerian one f2/f , 

n = n K =(^fj'\ (i) 

Note, just before a last stable orbit around a black hole, and of course inside 
it, this suggestion fails, but in the "standard" accretion disc model the relation 
(Q) is suggested to be fulfilled all over the disc, with an inner boundary at the 
last stable orbit. 

The equilibrium equation in the vertical z-direction is determined by a bal- 
ance between the gravitational force and pressure gradient 

dP GMz 

= (2) 

For a thin disc this differential equation is substituted by an algebraic one, 
determining the half-thickness of the disc in the form 

. i / p\ 1/2 



The balance of angular momentum, related to the <f> component of the Euler 
equation has an integral in a stationary case written as 

M{j - j in ) = -2nr 2 2ht rt p, t T<t> = (4) 

Here j = v^r = fir 2 is a specific angular momentum, t r( f, is a component of the 
viscous stress tensor, M > is a mass flux per unit time into a black hole, jo is 
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an integration constant having, after multiplication by M, a physical scncc of 
difference between viscous and advective flux of the angular momentum, when 
ji n itself is equal to the specific angular momentum of matter falling into a black 
hole. In the standard theory the value of ji„ is determined separately, from 
physical considerations. For the accretion into a black hole it is suggested, that 
on the last stable orbit the gradient of the angular velocity is zero, corresponding 
to zero viscous momentum flux. In that case 

j in = VL K r 2 ln , (5) 

corresponding to the Keplerian angular momentum of the matter on the last 
stable orbit. During accretion into a slowly rotating star which angular velocity 
is smaller than a Keplerian velocity on the inner edge of the disc, there is a 
maximum of the angular velocity close to its surface, where viscous flux is zero, 
and there is a boundary layer between this point and stellar surface. In that 
case (||) remains to be valid. The situation is different for accretion discs around 
rapidly rotating stars with a critical Keplerian speed on the equator. Here there 
is no extremum of the agular velocity of the disc which smoothly joins the star. 
In stationary self-consistent situation when the accreting star remains to rotate 
critically during the process of a disc accretion, the specific angular momentum 
of matter joining the star is determined by a relation (Bisnovatyi-Kogan, 1993): 
ji n = jjjj\crit: where the derivative is taken along the states of the star having a 
Keplerian equatorial speed. For stars with a polytropic structure, corresponding 
to equation of state P = Kp 1+ ~ , this derivative is calculated numerically giving 
the value ji n = 0.176fi,K-rf„ for n — 1.5; for n — 2.5; and negative values of 
jin for larger n. 

Note, that in the pioneering paper of Shakura (1972) the integration con- 
stant ji n was found as in (|^), but was taken zero in his subsequent formulae. 
Importance of using ji n in the form (||) was noticed by Novikov and Thorne 
(1973), and became a feature of the standard model. 

2.2 Viscosity 

The choice of the viscosity coefficient is the most difficult and speculative prob- 
lem of the accterion disc theory. In the laminar case of microscopic (atomic or 
plasma) viscosity, which is very low, the stationary accretion disc must be very 
massive and very thick, and before its formation the matter is collected by disc 
leading to a small flux inside. It contradicts to observations of X-ray binaries, 
where a considerable matter flux along the accretion disc may be explained only 
when viscosity coefficient is much larger then the microscopic one. In the paper 
of Shakura (1972) it was suggested, that matter in the disc is turbulent, what 
determines a turbulent viscous stress tensor, parametrized by a pressure 

t rr f, = -apv 2 = -aP, (6) 
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where v s is a sound speed in the matter. This simple presentation comes out 
from a relation for a turbulent viscosity coefficient rjt ~ pvtl with an average 
turbulent velocity vt and mean free path of the turbulent element /. It follows 
from the definition of t rr j, in ([!]), when we take I w h from (^) 

c2Q 2 . . 

t r< j, = pv t hr— fts pv t v a = —apv s , (7) 

where a coefficient a < 1 is connecting the turbulent and sound speeds v t = av s . 
Presentations of t rc f, in (g) and (0) are equivalent, and only when the angular 
velocity differs considerably from the Keplerian one the first relation to the right 
in (0) is more preferable. That does not appear (by definition) in the standard 
theory, but may happen when advective terms are included. 

Developement of a turbulence in the acctretion disc cannot be justified sim- 
ply, because a Keplerian disc is stable in linear approximation to the develope- 
ment of perturbations. It was suggested by Ya.B.Zeldovich, that in presence of 
very large Reynolds number Re = ^- the amplitude of perturbations at which 
nonlinear effects become important is very low, so in this situation a turbulence 
may develope due to nonlinear instability even when the disc is stable in linear 
approximation. Another source of viscous stresses may arise from a magnetic 
field, but it was suggested by Shakura (1972), that magnetic stresses cannot 
exceed the turbulent ones. 

Magnetic plasma instability as a source of the turbulence in the accretion 
discs has been studied extensively in last years (see review of Balbus and Hawley, 
1998). They used an instability of the uniform magnetic field parallel to the axis 
in differentially rotating disc, discovered by Velikhov (1959). It could be really 
important in absence of any other source of the turbulence, but it is hard to 
belive that there is no radial or azimuthal component of the magnetic field in 
matter flowing into the accretion disc from the companion star. In that case the 
field amplification due to twisting by a differential rotation take place without 
necessity of any kind of instability. 

It was shown by Bisnovatyi-Kogan and Blinnikov (1976, 1977), that inner 
regions of a highly luminous accretion discs where pressure is dominated by 
radiation, are unstable to vertical convection. Developement of this convection 
produce a turbulence, needed for a high viscosity. Other regions of a standard 
accretion disc should be stable to developement of a vertical convection, so other 
ways of a turbulence exitation are needed there. With alpha- prescription of 
viscosity the equation of angular momentum conservation is written in the plane 
of the disc as 

M(j-j in )=iTrr 2 aP h. (8) 

When angular velocity is far from Keplerian the relation (0) is valid with a 
coefficient of a turbulent viscosity 

V = apov s oh, (9) 
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where values with the index "0" denote the plane of the disc. 



= ^r— = -M- r [l-^). (10) 



2.3 Heat balance 

In the standard theory a heat balance is local, what means that all heat produced 
by viscosity in the ring between r and r + dr is radiated through the sides of 
disc at the same r. The heat production rate Q+ related to the surface unit of 
the disc is written as 

dVL 3 -GM ( %, 
'— = — M — g- 1 - — 
ar 87r r \ j 

Heat losses by a disc depend on its optical depth. The first standard disc model 
of Shakura (1972) considered a geometrically thin disc as an optically thick 
in a vertical direction. That implies enegry losses Q- from the disc due to a 
radiative conductivity, after a substitution of the differential eqiation of a heat 
transfer by an algebraic relation 

4 ncT 4 

Here a is a constant of a radiation energy density, c is a speed of light, T is 
a temperature in the disc plane, k is a matter opacity, and a surface density 
£ = 2ph. Here and below p, T, P without the index "0" are related to the disc 
plane. The heat balance equation is represented by a relation 



Q+ = Q-, (12) 

A continuity equation in the standard model of the stationary accretion flow is 
used for finding of a radial velocity v r 

v r = — = . (13) 

Anrhp 27rrE 

Equations (|l|),@,@, (p"2|), completed by an equation of state P(p, T) and rela- 
tion for the opacity k = n(p, T) represent a full set of equatiions for a standard 
disc model. For power low equations of state of an ideal gas P = P g = pTZT 
{JZ is a gas constant), or radiation pressure P — P r = , and opacity in 
the form of electron scattering n e , or Karammers formulae Kk , the solution of a 
standard disc accretion theory is obtained analytically (Shakura, 1972; Novikov, 
Thorne, 1973; Shakura, Sunyaev, 1973). Checking the suggestion of a large op- 
tical thickness confirms a self-consistency of the model. One of the shortcoming 
of the analytical solutions of the standard model lay in the fact, that solutions 
for different regions of the disc with different equation of states and opacities 
are not matched to each other. 
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2.4 Optically thin solution 

Few years after appearence of the standard model it was found that in addition 
to the opically thick disc solution there is another branch of the solution for 
the disc structure with the same input parameters M, M, a which is also self- 
consistent and has a small optical thickness (Shapiro, Lightman, Eardley, 1976). 
Suggestion of the small optical thickness implies another equation of energy 
losses, determined by a volume emission Q_ « q ph, where due to the Kirghoff 
law the emissivity of the unit of a volume q is connected with a Plankian aver- 
aged opacity k p by an approximate relation q « acT^Kp. Note, that Krammers 
formulae for opacity are obtained after Rosseland averaging of the frequency 
dependent absorption coefficient. In the optically thin limit the pressure is de- 
termined by a gas P = P g . Analytical solutions are obtained here as well, from 
the same equations with volume losses and gas pressure. In the optically thin 
solution the thickness of the disc is larger then in the optically thick one, and 
density is lower. 

While heating by viscosity is determined mainly by heavy ions, and cooling 
is determined by electrons, the rate of the energy exchange between them is 
important for a thermal structure of the disc. The energy balance equations are 
written separately for ions and electrons. For small accretion rates and lower 
matter density the rate of energy exchange due to binary collisions is so slow, 
that in the thermal balance the ions are much hotter then the electrons. That 
also implies a high disc thickness and brings the standard accretion theory to 
the border of its applicability. Nevertheless, in the highly turbulent plasma 
the energy exchange between ions and electrons may be strongly enhanced due 
to presence of fluctuating electrical fields, where electrons and ions gain the 
same energy. In such conditions difference of temperatures between ions and 
electrons may be negligible. Regretfully, the theory of relaxation in the turbulent 
plasma is not completed, but there are indications to a large exhancement of 
the relaxation in presence of plasma turbulence, in comparison with the binary 
collisions (Quataert, 1997). 

2.5 Accretion disc structure from equations describing 
continuously optically thin and optically thick disc re- 
gions 

In order to find equations of the disc structure valid in both limiting cases of 
optically thick and optically thin disc, and smootly describing transition between 
them, Eddington approximation had been used for obtaining formulae for a 
heat flux and for a radiation pressure (Artemoma et al., 1996). The following 
expressions had been obtained for the vertical energy flux from the disc Eq, and 
the radiation pressure in the symmetry plane 
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2acT n 4 A 4 2 \ 1 „ aT 4 1 + 3^ 



F o = ^"^ 1 + ^- + ^ , ^,0 = ^ - , 4 2 . (14) 
3r V 3r 3r2 / 3 1 + * + 2 



3to 3r^f 

where r = « e p/i, r* = (toTqo) 1 ^ 2 , t q0 « At r 3> t* 3> 1 we have (O) 

from (Jl4]). In the optically thin limit r, « t < 1 we get 

2 

F = acT r Q0 , -FWo = -acT Q T a Q. (15) 

Using Fb instead of Q- and equation of state P = pTZT + P ra d,o, the equations 
of accretion disc structure together with equation Q+ = Fo, with Q+ from (|l0|), 
have been solved numerically by Artemova et al. (1996). It occures that two 
solutions, optically thick and optically thin, exist separately when luminosity 
is not very large. Two solutions intersect at m = rhb and there is no global 
solution for accretion disc at rh > rrib (see Fig.l). It was concluded by Artemova 
et al (1996), that in order to obtain a global physically meaningful solution at 
m > rhf,, account of advection is needed. 



3 Accretion discs with advection 

Standard model gives somewhat nonphysical behaviour near the inner edge of 
the accretion disc around a black hole. For high mass fluxes when central re- 
gions are radiation-dominated (F w P r , k w K e ), the radial dependence follows 
relations (Shakura, Sunyaev, 1973) 

p ~ r 3 / 2 J- 2 -> 00, T ~ r- 3 / 8 , (16) 
h~J^0, E ~ r 3/2 J- 1 ->• 00, u r ~ r~ 5/2 J -> 0, 
where limits relate to the inner edge of the disc with r = ri ni J = 1 — = 

1 — \/^F- At smaller M, when near the inner edge P w F g , k ps k g , there are 
different type of singularities 

p ~ r -33/20 j2 /5 ^ Qj T ^ r -9/10j2/5 ^ 0; 
ft „ _^ Q) j, _ r -3/B j3 /5 ^ Qj ^ r -2/5j-3/5 ^ ^ 

This results from the local form of the equation of the thermal balance (|l2| ) . It 
is clear from physical ground, that when a local heat production due to viscosity 
goes to zero, the heat brought by radial motion of matter along the accretion 
disc becomes more important. In presence of this advective heating (or cooling 
term, depending on the radial entropy S gradient) written as 

M dS 

Qadv = ~ T—, (18) 

2-Kr dr 
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Figure 1: The dependences of the optical depth To on radius, r* = r/r g , for 
the case Mbh = 10 8 M Q , a = 1.0 and different values of m. The thin solid, 
dot-triple dash, long dashed, heavy solid, short dashed, dotted and dot-dashed 
curves correspond to m = 1.0,3.0,8.0,9.35,10.0,11.0,15.0, respectively. The 
upper curves correspond to the optically thick family, lower curves correspond 
to the optically thin family. 
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the equation of a heat balance is modified to + Q a dv = Q-- In order 
to describe self-consistently the structure of the accertion disc we should also 
modify the radial disc equilibrium, including pressure and inertia terms 

r(fl 2 -n 2 K ) = -— -v r -^. (19) 
p dr dr 

Appearence of inertia term leads to transonic radial flow with a singular point. 
Conditions of a continious passing of the solution through a critical point choose 
a unique value of the integration constant ji n . First approximate solution for 
the advective disc structure have been obtained by Paczynski and Bisnovatyi- 
Kogan (1981), but a corresponding set of equations had beed discussed earlier 
(Hoshi and Shibazaki, 1977; Liang and Thompson, 1980). Attempts to find a so- 
lution for advective disc structure (see e.g. Matsumoto et al., 1984; Abramovicz 
et al., 1988) gave the following results. For moderate values of M a unique 
continuous transonic solution was found, passing through singular points, and 
corresponding to a unique value of ji n . The number of critical ponts in the 
radial flow happens always to be more then unity. This is connected with two 
reasons. First, the gravitational potential <p g in papers dealing with advective 
disc solutions was different from newtonian one (Pacziriski and Wiita, 1980): 
4>g = ~ i r g — 2G ^ L ■ The advantage of this potential is a realistic approx- 
imation of the general relativistic (GR) effects, namely, infinitive gravitational 
attraction at a gravitational radius r g , and existence of the stable circular orbits 
only up to r = 3r 9 , like in exact GR. Appearence of two critical points for a 
radial flow in this potential was analysed by Chakrabarti and Molteni (1993). 
The second reason of multiplicity of singular points is connected with using of 
equations averaged over a thickness of the disc. That changes a structure of 
hydrodynamic equations, leading to a position of singular points not coinsiding 
with a unit Mach number point, and increasing a number of critical points. 

When M is becoming so high, that radiation pressure starts to be important, 
still unresolved problems appear in a construction of the advective disc model. 
These problems are connected with increasing of a number of a critical points 
from one side, and loss of uniqueness of the transonic solution from another. 
So, with increasing of M the solution becomes nonunique at some parameters, 
or was not found at all (see Matsumoto et al., 1984; Abramovicz et al., 1988; 
Artemova et al., 1996a). At high M the integral curves are very sensitive to 
input conditions: form of viscosity stresses (0) or @, choice of boundary con- 
ditions etc. The system of equations has a very small resource of stability, so 
it cannot be excluded, that the failures are connected with an inproper choice 
of a numerical method and developement of numerical instabilities prevents of 
finding a unique physical solution. In addition to continuous solutions, solutions 
with standing shock waves have been investigated (Chakrabarti, 1996). 
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3.1 Two-temperature advective discs 



In the optically thin accretion discs at low mass fluxes the density of the matter 
is low and energy exchange between electrons and ions due to binary collisions 
is slow. In this situation, due to different mechanisms of heating and cooling for 
electrons and ions, they may have different temperatures. First it was realized 
by Shapiro, Lightman, Eardley (1976) where advection was not included. It 
was noticed by Narayan and Yu (1995), that advection in this case is becoming 
extremely important. It may carry the main energy flux into a black hole, 
leaving rather low efficiency of the accretion up to 1CP 3 — 1CP 4 (advective 
dominated accretion flows - ADAF). This conclusion is valid only when the 
effects, connected with magnetic field annihilation and heating of matter due to 
it are neglected. 

In the ADAF solution the ion temperature is about a virial one kTi ~ 
GMrrii/r, what means that even at high initial angular momentum the disc 
becomes very thick, forming practically a quasi-spherical accretion flow. It is 
connected also with an " alpha" prescription of viscosity. At high ion tempera- 
tures, connected with a strong viscous heating, the ionic pressure becomes high, 
making the viscosity very effective. So, due to suggestion of "alpha" viscos- 
ity in the situation, when energy losses by ions are very low, some kind of a 
"thermo- viscous" instability is developed, because heating increases a viscos- 
ity, and viscosity increases a heating. Developement of this instability leads to 
formation of ADAF. 

A full account of the processes, connected with a presence of magnetic field 
in the flow, is changing considerably the picture of ADAF. It was shown by 
Schwarzman (1971), that radial component of the magnetic field increses so 
rapidly in the spherical flow, that cquipartition between magnetic and kinetic 
energy is reached in the flow far from the black hole horizon. In the region 
where the main enegry prodiction takes place, the condition of cquipartition 
takes place. In presence of a high magnetic field the efficiecy of a radiation 
during accretion of an interstellar matter into a black hole increase enormously 
from <~ 1CP 8 up to <~ 0.1 (Schwarzman, 1971), due to efficiency of a magneto- 
brcmstrahlung radiation. So possibility of ADAF regime for a spherical accre- 
tion was noticed long time ago. To support the condition of equipartition a 
continuous magnetic field reconnection is necessary, leading to annihilation of 
the magnetic flux and heating of matter due to Ohmic heating. It was obtained 
by Bisnovatyi-Kogan and Ruzmaikin (1974), that due to Ohmic heating the 
efficiency of a radial accretion into a black hole may become as high as ~ 30%. 
The rate of the Ohmic heating in the condition of cquipartition was obtained 
in the form 



In the supersonic flow of the radial acrretion equipartition between magnetic 




3 B 2 



(20) 
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and kinetic energy was suggested by Schwarzman (1971): 



&_ „ pvl = pGM 

8tt ~ 2 r [ ' 

For the disc accretion, where there is more time for a field dissipation, almost 
cquipartition was suggested (Shakura, 1972) between magnetic and turbulent 
energy, what reduces with account of "alpha" prescription of viscosity to a 
relation 

8^ ~ — = 2 a ™ P ' (22) 
where a m characterises a magnetic viscosity in a way similar to the turbulent a 
viscosity. It was suggested by Bisnovatyi-Kogan and Ruzmaikin (1976) the sim- 
ilarity between viscous and magnetic Reynolds numbers, or between turbulent 
and magnetic viscosity coefficients 

pvl pvl 
Re = , Re m = , (23) 

V Vm 

where the turbulent magnetic viscosity rj m is connected with a turbulent con- 

2 

ductivity a — ^ . Taking r\ m — ^r], we get a turbulent conductivity 

in the optically thin discs. For the radial accretion the turbulent conductivity 
may contain mean free path of a turbulent element l t , and turbulent viscocity 
v t in (HJ) instead of h and v s . In ADAF solutions, where ionic temperature is 
of the order of the virial one the two suggestions (|2l|) and (^2|) almost coinside 
at a m — 1. 

The heating of the matter due to an Ohmic dissipation may be obtained 
from the Ohm's law for a radial accretion in the form 

dS a£* ^ vlB* B*v t 

dr v r v r c^ / vKa m v r lt 

what coinsides with (|20|) when a m = Jp^- , or l t = ^ rv ^ . Here a local electrical 
field strength in a highly conducting plasma is of the order of £ ~ for the 
radial accretion. 

Equations for a radial temperature dependence in the accretion disc, separate 
for the ions and electrons are written as (Bisnovatyi-Kogan and Lovelace, 1997) 

^ = 'H-qi + Ti-Bi — Qie , (26) 



dt p 
dE e P e dp 
~dT ~ ~~{?~dt 



f^rje ~t~ ^Se ~t~ Qie ^brem &cyc ? 

(27) 
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Here 4t = St Vr Sr' ^ ra ^ e °^ a v i scous heating of ions 7i ri i is obtained from 
© as 

-w 27rr n 3 ~ VKV s n-i [™e<y (or\ 

M 2 r V m i 

where vk = tSIk- The rate of the energy exchange between ions and electrons 
due to the binary collisions was obtained by Landau (1937), Spitzer (1940) as 



4(27r)2ne 4 ( T e T, 



m,im e \m e m 



inA(Ti-T e ), (29) 



with £nA = O(20) the Coulomb logarithm. The electron bremstrahlung Cbrem 
and magneto-bremstrahlung C cyc cooling are taken into account. The expression 
for an Ohmic heating in the turbulent accretion disc may be written in different 
ways, using different velocities ve in the expression for an effective electrical 
field £ = 2Ls5, A self-consistency of the model requires, that expressions for a 
magnetic heating of the matter Hb > obtained from the condition of stationarity 



of the flow (p0|), and from the Ohm's law (25), should be identical. That gives 



some restrictions for the choice of a characteristic velocity ve- Comparison 



between (20) and ( |25|) shows the identity of these two expressions at ve = 
v r , , = So, the model is becoming self-consistent at the reasonable 
choice of the parameters. Note, that in the advective models J is substituted by 
another function which is not zero at the inner edge of the disc. The heating due 
to magnetic field reconnection Hb in the equations (p6|), (p7|), may be written 
as 

167T rp 2 J \vk J 

So, at vb = vk the expressions for viscous and magnetic heating are almost 
identical. The distribution of the magnetic heating between electrons and ions 
has a critical influence on the model, if we neglect the influence of a plasma 
turbulence on the energy relaxation, and take into account only the energy 
exchange by binary collisions from (p9|). Observations of the magnetic field 
reconnection in the solar flares show (Tsuneta, 1996), that electronic heating 
prevails. 

It follows from the physical picture of the field reconnection, that transfor- 
mation of the magnetic energy into a heat is connected with the change of the 
magnetic flux, generation of the vortex electrical field, accelerating the parti- 
cles. This vortex field has a scale of the turbulent element and suffers rapid 
and chaotic changes. The accellerating forces on electrons and protons in this 
fields are identical, but accelerations themselfs differ ~ 2000 times, so during 
a sufficiently short time of the tubulcnt pulsation the electron may gain much 
larger energy, then the protons. Additional particle acceleration and heating 
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happens on the shock fronts, appearing around turbulent cells, where reconnec- 
tion happens. In this process acceleration of the electrons is also more effective 
than of the protons. In the paper of Bisnovatyi-Kogan and Lovelace (1997) the 
equations (^), ( p7| ) have been solved in the approximation of nonrelativistic 
electrons, vb=vk, what permitted to unite a viscous and magnetic heating into 
a unique formila. The combined heating of the electrons and ions were taken 
as H e = (2 — g)TC V i, TC e = gTC V i- In the expression for a cyclotron emission 
self-absorption was taken into account according to Trubnikov (1973). The re- 
sults of calculations for g = 0.5 1 show that almost all energy of the electrons 
is radiated, so the relative efficiency of the two-temperature, optically thin disc 
accretion cannot become lower then 0.25. Note again that accurate account of 
a plasma turbulence for a thermal relaxation and corresponding increase of the 
term Qi e may restore the relative efficiency to its unity value, corresponding to 
the optically thick discs. 

4 Discussion 

Observational evidences for existence of black holes inside our Galaxy and in the 
active galactic nuclei (Cherepashchuk, 1996, Ho, 1998) make necessary to revise 
theoretical models of the disc accretion. Large part of high energy radiation 
indicates to its origin close to the black hole, where standard accretion disc 
model is not a appropriate. The improvements of a model are connected with 
account of advective terms and more accurate treatement of the magnetic field 
effects. Conclusions about existence of ADAF solution for an optically thin 
accretion disc at low mass flux are connected with an incomplete account of 
the effects connected with magnetic field annihilation. Their account does not 
permit to make a relative efficiency of the accretion lower then ~ 0.25 from the 
standard value. It is expected that more accurate treatement of the relaxation 
connected with the plasma turbulence will even more increase the efficiency, 
making it close to unity (see also Fabian and Rees, 1995). 

Some observational data which were interpreted as an evidence for the ex- 
istence of the ADAF regime have disappered after additional accumulation of 
data. The most interesting example of this sort is connected with the claim 
of the proof of the existence of event horizon of the black holes due to mani- 
festation of the ADAF regime of accretion (Narayan et al., 1997). Analysis of 
the more complete set of the observational data (Chen et al., 1997) had shown 
disapperence of the statistical effect clamed as an evidence for ADAF. This ex- 
ample shows how dangerous is to base a proof of the theoretical model on the 
preliminary observational data. It is even more dangerous, when the model is 
physically not fully consistent. Then even a reliable set of the observational data 
cannot serve as a proof of the model. The classical example from astrophysics of 
this kind gives the theory of the origin of the elements presented in the famous 
book of G.Gamov (1952), where the model of the hot universe was developed. In 
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addition to rich advantages of this model, the author also wanted to explained 
the origin of heavy elements in the primodial explosion, neglecting the problems 
connected with an absence of the stable elements with the number of barions 
equal to 5 and 8. G. Gamov considered a good coinsidence of his calculations, 
where the mentioned problem was neglected, and the observational curve, as a 
best proof of his theory of the origin of the elements. The farther developements 
have shown that his outstanding theory explains lot of things, exept the origin 
of the heavy elements, which are produced due to stellar evolution. 

It looks like it is difficult to use ADAF for solution of the problem of existence 
of undcrluminous AGN, where the observed flux of the energy is smaller, then 
the expected from the standard accretion disc models. Two possible ways may 
be suggested. One is based on a more accurate estimations of the accretion 
mass flow into the black hole, which could be overestimated. Another, more 
attractive possibility, is based on existence of another mechanisms of the energy 
losses in the form of accelerated particles, like in the radio-pulsars, where their 
losses exceed strongly a radiation losses. This is very probable to happen in 
a presence of a large scale magnetic field which may be also responsible for 
a fomation of the observed jets. To extend this line, we may suggest, that 
underlumilnous AGN loose main part of their enegry to the formation of jets. 
The search of the correlation between existence of jets and lack of the luminosity 
could be very informative. 
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